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Abstract. In this paper we evaluate sources of background for the 7p — ► U)p, with the u) detected through its 
7r°7 decay channel, to compare with the experiment carried out at ELSA. We find background from 7p — > 
TT°ir°p followed by decay of a 7r° into two 7, recombining one 7r° and one 7, and from the 7p —¥ iv rip reaction 
with subsequent decay of the rj into two photons. This background accounts for the data at 7r°7 invariant 
masses beyond 700 MeV, but strength is missing at lower invariant masses which was attributed to photon 
misidentification events, which we simulate to get a good reproduction of the experimental background. 
Once this is done, we perform an event mixing simulation to reproduce the calculated background and we 
find that the method provides a good description of the background at low 7r°7 invariant masses but fakes 
the background at high invariant masses, making background events at low invariant masses, which are 
due to 7 misidentification events, responsible for the background at high invariant masses which is due to 
the 7p — > ty°ty°p and 7p — > -K°rip reactions. 

PACS. 13.60.Le describing text of that key - 25.20.Lj describing text of that key 



1 Introduction 

The interaction of vector mesons with nuclear matter has 
attracted attention for a long time and has been tied to 
fundamental aspects of QCD. Yet, the theoretical mod- 
els offer a large variety of results from a large attraction 
to a large repulsion. Early results on this issue within 
the Nambu Jona Lasinio model produced no shift of the 
masses pQ while, using qualitative arguments, a universal 
large shift of the mass was suggested in [2]. More recent 
detailed calculations show no shift of the mass of the p 
meson in matter [3U31[3] . Experimentally the situation has 
undergone big steps recently with the NA60 collaboration 
reporting a null shift of the p mass in the medium [6,7 
in the dilepton spectra of heavy ion reactions and also a 
null shift in the 7 induced dilepton production at CLASS 
[5] . On the other hand the KEK team had earlier reported 
an attractive mass shift of the p in [51HU]. As explained 
in detail in [5], the different conclusions can be traced 
back to the way the background is subtracted. Thus, the 
treatment of the background is an essential part of the in- 
vestigation of the vector meson properties in nuclei. The 
case of the uj in the medium is more obscure. Theoreti- 
cally there are about twenty different works with claims 
from large attraction to large repulsion (see [PTlll"2"llT3] for 
details). Experimentally there are claims of a large shift 



of the mass of the uj from the study of the photon induced 
uj production in nuclei, with the uj detected through its 
7r°7 decay channel [141IT5] . However, it was shown in [T^] 
that the shift could just be a consequence of a particu- 
lar choice of background subtraction and that other rea- 
sonable choices led to different conclusions. For instance, 
choosing a background in the nuclear case proportional to 
the background on the proton in the region below the uj 
peak, the experimental data could be explained without a 
shift of the u) mass in the medium. 

The method to determine the mass shift is very dif- 
ferent from the one used to determine the width in the 
medium. This latter one relies upon the production cross 
section in different nuclei, which leads to the transparency 
ratio that allows to determine a large width of the uj in 
the medium |12U16j . By contrast the measurement of the 
mass requires the analysis of the shape of the invariant 
mass distribution, which is barely affected in nuclei be- 
cause practically all u) decays occur in the nuclear surface 
or outside the nucleus. 

The discussion on the issue of the uj mass shift was 
followed by the evaluation of the background with the 
mixed event technique in |17j . There the background for 
a nuclear target was evaluated and found to be the same 
as assumed in [TS], and again it was concluded that the 
data demanded a shift of the uj mass in the medium. 
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The former discussion indicates that the treatment of 
background is an essential issue in this problem. In view 
of this we decided to face the problem and investigate the 
details on how the mixed event technique works in the 
present case. For this we followed the strategy of evaluat- 
ing the background for the proton target. We could trace 
two sources of background that account for the experi- 
mental cross section at 7r°7 invariant masses of the order 
of the uj mass and beyond. The rest of the background 
at lower invariant masses was simulated to account for 7 
misidentification events, as found in |14) . Once a back- 
ground consistent with the experiment is obtained theo- 
retically, we apply the mixed event technique to obtain 
the background and compare it with the theoretical one. 
In this way we can determine the ability of the mixed event 
method to reproduce the background in this reaction. The 
results that we find are that the method can provide the 
background at low invariant masses, where the cross sec- 
tions are large, but it actually fakes the background in the 
region of invariant masses around and beyond the uj mass. 
We show that the mixed event generated background in 
that region is completely tied to the real events at low 
7r°7 invariant masses where the da / dM^o 1 distribution is 
larger. As a result, we show that the distribution obtained 
with the mixed event method in the region of invariant 
masses around and beyond the uj mass is largely insensi- 
tive to the actual background contributing in that region. 
But this is precisely the region where the background is 
needed to determine changes of the uj signal in the nu- 
clear medium. As a consequence we clearly show that the 
mixed event technique is unsuitcd in the present case as 
an instrument to determine possible shifts of the uj mass 
in the medium. 

In the meantime, a recent reanalysis of the background 
of the reaction of [14irT5] done in [18] concludes, however, 
that one cannot claim a shift of the uj mass from this 
experiment. 



2 Background sources in the reaction 

7P — > up (ui — > 7T°7) 

The reaction that we study is jp — > ujp where the uj is 
detected through its uj — > 7r°7 decay mode. This is the 
reaction studied in the CBELS A/TAPS experiment. Ac- 
cording to the study in [14], in the region of the recon- 
structed invariant mass of the 7r°7 one of the main sources 
of background comes from the jp — » ir°ir°p reaction fol- 
lowed by the decay of any of the two ir° into 77. Then, 
background events appear from the combination of one 
of these photons and the remaining ir° . Another source is 
the 7p — !• TT°rjp reaction. We evaluate the cross section for 
those two processes in the following subsection. 



2.1 The 7p — > ir°ir°p reaction 

This reaction has been thoroughly studied at Mainz 20, 
|2T J 2"2 1| 2"5] and more recently at EISA [HES], GRAAL [25] 



and Jefferson Lab [2Tl[281f25] . We are interested not only in 
the cross section for the reaction but at the same time to 
have an event generator that provides events weighed by 
their probability determined by the available phase space. 
For this purpose the Monte Carlo evaluation of the cross 
section integral is the most suited algorithm since it pro- 
vides the events allowed by phase space properly weighted 
and the cross section in the end. 

The 7p — > 7r 7r°p cross section is given by 

= M 2 N f d 3 Pl f d 3 P2 f d 3 p 3 
a s-M 2 N J (2tt)3 J (2tt) 3 J (2tt)3 



2E{ Pl ) 2E{p 2 ) 2E{p 3 ) 1 1 
x(2tt) 4 60 7 +p p -pi -p 2 -p 3 ), (1) 

which includes the 1/2 symmetry factor to account for 
the two identical pions in the final state. The variables 
Pi, P2, P3 are the momenta of the final proton and the two 
7T° respectively and T stands for the transition matrix for 
the 7p — > tt q tt°p process. In the \T\ 2 factor a sum over final 
spins and a proper average over initial spins is implicit. 

Since d 3 p/E(p) is a Lorentz invariant measure, we pro- 
ceed to evaluate the last two integrals in Eq. ([1]) in the 
reference frame where p 7 + p p — p\ = which guarantees 
that P2 + p 3 = 0. The cross section is then written as 

where M23 is the invariant mass of the two 7r° given by 
M2 3 = (p~ ( +P P ~Pi) 2 = s + M 2 N - 2{ Pl +p p )p 1 . (3) 

We shall work in the laboratory (lab) frame that allows us 
to implement easily all the experimental cuts. The variable 
p 2 in Eq. ^ is the ir° momentum in the 7r°7r° rest frame 

- A^(M| 3 ,m^,m^) 

P2 2M2-3 ' W 

and dfl is performed in the 7r°7r° rest frame. 

The jp — > ir°ir°p reaction has been studied theoret- 
ically at E 1 < 0.8 GeV J301I3I1E21E31E1] and at higher 
energies E 1 > 1.5 GeV in [27,35]- We do not need any of 
these sophisticated models here. Our final goal is to see 
how much background comes form this source and to have 
an event generator by means of which we can study how 
the mixed event technique works in the present case. For 
this purpose it is enough to consider |T| 2 to be a constant 
over the phase space and fit its value to reproduce the 
experimental results for the cross section in the region of 
interest to us. We take the cross section for ^p — > 7r°7r°p 
at the needed photon energies from [24lf25] . 

The next step is to write p~2 in the lab frame. We have 
in the 7r°7r° rest frame 

!sin9cosip J 
sin9 sirup > (5) 
cos9 
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P3 



P2 



(6) 



with 9, ip angles in the 7r°7r° rest frame. We then perform a 
boost of f>2 to the lab frame where p 2 +p 3 = p 7 +p p Pi = 



P2 



( E 2 



2.; 



P2 



P°2 



p2 



M 



2;s 



P+P2, 



(7) 



where E 23 is the two pion energy in the lab frame E 23 — 
(Af 2 2 3 + P 2 )3. Similarly we boost p 3 to p 3 in the initial jp 
lab frame. 

Assume now that the pion with momentum p 2 is the 
one that decays into 77. In the pion rest frame the two 7's 
will go back to back and one 7 will have the momentum 



p 7 



sm9~ l cosip-y 
sin9 1 sirnp 1 
cos9^ 



(8) 



with 8 ' 7 , ip j angles of the photon in this one ir° rest frame. 
Once again we boost this photon momentum to the frame 
where the pion has momentum p 2 



p 7 



E-2 



- 1 



Pi P2 
Pi 



P2+P 



7' 



0) 



Since we can have a 7r°7 combination from either of the 
two 7T° or the two 7's, we would obtain a combinatorial 
factor of four to account for these possibilities. 

All this said, the cross section for jp — > tt q tt°p —> 
77 ir°p — > 77T° + X reads 



<7 = 4 



71/ 



N 



1 \T\ 2 9(M 23 - 2m,o) (10) 



M%J (2^)3 2%) 



1 ,2tt 111/" 1 

< I dcorij Apfr——-] 



2tt 

dcos9 1 I dp> 1 



recalling from Eq. ([3]) that in the lab frame 

M 23 = s + M% - 2{E lin + M N )E{ Pl ) + 2p jin ■ pi(ll) 

Next one generates random numbers for pi , 9, ip, 6* 7 , ip 7 
with \pi I restricted between zero and 



\Pl\lab 



|„ \max 1 Tpmax 



(12) 



with and i?] 7 ^^ the maximum momentum and 

energy allowed for the final proton in the jp center of mass 
(CM) frame, corresponding to the case when the two 7r° 
go together 



\Pl\CM - — 



(13) 



and v is the velocity of the jp CM system measured in 
the lab frame 



\Pi 



E 7 in + Mpf 



(14) 



For each of these events we evaluate the invariant mass 
M int ,(7r°7) = (p 7 +p 3 ) 2 (15) 

and store the events, properly weighted by |T| 2 and phase 
space factors, in boxes of Mi nv (ir°j) for a suitable parti- 
tion of M inv (w j). 



2.2 The 7P — > 7r°?7p reaction 

We also evaluate the background for the jp — > ■K A r\p reac- 
tion. In this case r\ — >■ 77 and we get one photon from there 
plus the 7T° to reconstruct the 7r°7 invariant mass. The 
changes with respect to the former reaction are minimal. 
Since now there is only one ir° , we do not have to include 
the 1/2 symmetry factor and the combinatorial factor of 
four before is now a factor of two and the mass of one pion 
must be changed to the mass of the rj when needed. There 
are data for this reaction in [36 , 37 ,38 ,39]. There are also 
recent detailed models for the reaction [4T)PT] accounting 
fairly well for the cross section |3"6"ll3"T] and the asymme- 
tries [38 . However, once again, for the present problem it 
suffices to repeat the procedure done for the 7^ — > 7r°7r°p 
reaction in the former section taking a constant |T| 2 and 
implementing properly the phase space demanding that 
we reproduce the data of [3B1l37U3tJj. 



2.3 Background from extra sources 

As we shall see in the results section, the 7^ — > ir°ir°p 
and 7p — > ir°r)p reactions can account for the background 
observed in the CBELSA/TAPS experiment [T4lfT5] in the 
region of the to excitation and higher 7r°7 invariant masses. 
However, it docs not account for the large background 
observed in the region of 7r°7 invariant masses lower than 
m^. Once again we resort to the findings of [T3] suggesting 
that such events could come from reactions like *yp — > 
7T 7r + n with a misidentification of the neutron by a photon 
and other possible sources of 7 misidentification. In this 
part we do not make a theory since the events come from 
ignorance of the occurring reactions and misidentification 
of particles which have to do with the detector system. 
However, we would like to have events corresponding 
to this region in order to perform later on the mixed event 
analysis. 

To generate background events in this region we write 
the corresponding cross section as 



dM ff7 |T(M W7 )| S 



(16) 



I dcos9 
-1 Jo 



2tt 



d V / d 3 P CM 6>(500McV/c 



CM 



I), 



where |T(M 7r7 )| 2 is a function to be determined from ex- 
periment. 8, ip are the tt° angles in the 7r°7 rest frame. 
There the pion momentum is 



sin9 cosip 
p^o — p^o ^ sin9 simp 
cos9 



(17) 
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with 



Pir« 



M 2 



2AU 



Pi 



-PlT<> ■ 



(18) 



Eq. (I16[) contains an integral over Pcm with a maximum 
of 500 MeV for I-PcmI- This momentum represents the 
7r°7 total momentum in the "fi n p CM frame. In this frame 
the momenta are more evenly distributed than in the lab 
frame and then we take an isotropic distribution for Pq m 
with the constraint |Pca/| < 500MeV/c. This is a conser- 
vative estimate that exceeds the phase space of Eq. JUT 



when \P\ = |jvi+p 7 | the total it 7 momentum in the "fi n p 
lab frame is restricted to values smaller than 500 MeV/c. 
Next we boost the ir° and the 7 momenta from their CM 
frame to the ji n p CM frame where the 7r°7 system has 
momentum Pcm- We have 



Pit" 



Ptt° 



CM 



rP 

Ptt'> 



p2 

CM 



Mir 



CM + Ptt« , 



^cm 



E 

pL = Pcm - Piro 



M 2 



(19) 
(20) 
(21) 



The next step is the boost to the lab system where the 
"finP has momentum p 7 i n and energy E lp — p lin + M N , 
hence 



Pir° 



E 

J-/>yp 

7^ 



P,o ■ Pi 



r>'° 



P 2 lln V~ s 



P-yin+PnO, (22) 



and a similar one for p 7 . On these 7 and 7r° momentum 
we enforce now the cut 



is particularly simple to implement. The Monte Carlo inte- 
grals are done by generating random events into a volume 
V containing the whole phase space and then the integral 
is given by the average value of the integrand, that we 
shall denote as |T| 2 , times the volume 



Eli \n 2 

N 



V. 



(25) 



where one understand that \Ti\ 2 is zero if the event gen- 
erated does not belong to the phase space. Assume now 
that we take pair of events i, j corresponding to a reaction 
channel. We have 



E 



mi 



N 



V 



E 



N 



i\Tj\ 2 



or equivalently 



N N 
i=\ j=l 



m 



N 



W\T 3 \ 2 



V, 



V 



N 2 



(26) 



(27) 



We can then generate pairs of events in the phase space 
volume V and the former integral gives the cross sec- 
tion. Simultaneously with the evaluation of the cross sec- 
tion we can obtain M w o 7 for each pairs of events as done 
in Eq. (|24j) and store the event, with its corresponding 
weight, in a box of a certain M^a^ value. After the dou- 
ble sum in Eq. (|2T|) we obtain the normalized da/dM n o^ 
distribution. 

The generalization to four channels as we have in our 
case, 7p — » 7r°7r°p, jp — » ir°r]p, the channel from 7 misiden- 
tification and the jp — > top — > ir°^p channel is straightfor- 
ward 



IPtt" +P 7 | < 500MeV/c. 



(23) 



The function T(M 7ri ) of Eq. (TTB1) is determined empirically 



n°w°p 



such that the sum of the cross section for jp 
plus 77) — » ir a r)p, plus the new one simulating 7 misidenti- 
fication events, gives the total experimental cross section 

of pins]. 



Ctot 



Otot 



N N 

EEEE 

a i—1 j—1 



where 



(l)l2 F | T jJ|2 
~2 



Otot — Ol + &2 + CT3 + (7 4 , 



V 



(28) 



(29) 



and a, (3 run from 1 to 4. In order to obtain da/dM^a 



3 Mixed events calculation 

In the mixed event simulation the idea is to obtain the 
background from the real data by evaluating M 1T a 1 se- 
lecting the 7T° and the 7 from two different events. The 
invariant mass distribution is then given by 



M 2 JME) = (p 7r o(l)+p J (2)) 2 



(24) 



There is abundant literature on the subject [421I43II44U451 
146) and it has become a popular instrument to determine 
background and isolate particular reactions that peak at 
a certain place. 

In our case where all the integrals of the cross sections 
are performed by Monte Carlo, the mixed event simulation 



from these mixed events we evaluate again M 7r o 1 



from 



Eq. (j2~4")l and store the events in boxes of M^o^ and we 
obtain at the end the histogram that provides us with the 
normalized da / dM^o 7 distribution from mixed events. 



4 Results 

In Fig.[T]the contributions from the cj-signal and different 
sources discussed above are compared with the experimen- 
tal invariant mass da / dM 7T a 1 distribution in the reaction 
7P 7r°7p from CBELSA/TAPS experiment. 

The sources of background are jp ~ > Tr°w°p with either 
of the two pions decaying into two 7, which was stud- 
ied in section 2.1. This is the most important source of 



5 




Fig. 1. The invariant mass do /dM K a~ distribution in the reaction 7p — > Tv°jp. The signal and different background contributions 
are also shown. 



background in the region of the lo and beyond. The other 
important source of background is the 7 misidentification 
studied in Section 2.3. This source competes with the for- 
mer one in the region of the omega and becomes dom- 
inant at smaller 771" invariant masses. The third source 
considered is the one coming from jp — > Tr°rjp followed 
by the r\ decay into two 7. This source has a smaller 
strength than the other two, but was found to be im- 
portant to understand a peak in the experiment at lower 
77T° invariant masses than the omega in [13] when protons 
were measured in coincidence. The omega signal comes 
from our study in [12]. The fit to the unnormalized data 
is done by adjusting the strength of the ■yp — > tt tt°p 
source to the experiment distribution at large invariant 
mass, since this is the most important source in this in- 
variant mass region. The source of the 7p — > ir°r]p, as well 
as the signal are rescaled keeping their ratio, in order to 
keep the theoretical proportion between all these differ- 
ent sources. Finally, the source of misidcntified 7 is added 
in order to complete a good description of the data. As 
one can see from the figure, the agreement of the theoret- 
ical model with the experimental data is very good. Note 
that we also have adapted our theoretical set up to the 
experimental one by choosing E 1 — 0.9 — 2.6 GeV and 
Ptt°7 = \Ptt° + P7I < 500 MeV. 



We should note that the experimental spectrum shown 
in Fig. [1] is not acceptance corrected. The reason is that 
experimentally one observes only three out of four photons 
in the final state due to the detector inefficiencies, or the 
overlap of photon clusters, and the latter depends on the 
energy of the photon [18]. Yet, in the unnormalized spec- 
trum to which we make the fit, the differences are of the 
order of 20 % from the lower mass part of the spectrum 
to the higher mass one, and have irrelevant consequences 
for the argumentations and conclusions that follow. The 
consideration of acceptance would be more important in 
case one would like to compare our different sources of 
background with the acceptance uncorrected experimen- 
tal determinations |18j . which is not our concern here. 



5 Mixed events and different tests 

On the first hand we make the ME simulation that was 
done in [17] taking two independent events and demanding 
that 

Ip^+p^l < 500MeV/c (30) 

for the mixed event. We call it method I. This choice 
has in principle a conceptual flaw. Indeed, the curve in 




M o [MeV] 

k y 

Fig. 2. The mixed event background with method I is shown, together with the actual background, the results of the model 
for ui signal plus background and the data. 



Fig. [TJ corresponds to events in which one has imposed 
Ip^o + p 7 | < 500MeV/c in the 7r° and 7 momenta of 
the same event. This restricts the phase space consider- 
ably. Now if one imposes Eq. after the mixing, one 
can have both events (1) and (2) or one of them that 
do not fulfill separately |p^o + p 7 | < 500 MeV/c and, as 
a consequence, these are events which do not contribute 
to the distribution of Fig. [TJ In other words, one can 
be using events that do not provide any information to 
Fig. [TJ to obtain its corresponding background through 
the mixed event method. Clearly in the extreme case that 
most events in the ME simulation do not pass the individ- 
ual \p w o + p 7 | < 500MeV/c test one would be obtaining 
the background of the curve from a physical situation that 
has no relationship with the distribution of Fig. [TJ How 
far is one in practice from this situation depends of course 
on the cut. 

In Fig. [5] we show the results that we get from the 
mixed event simulation for the background, compared to 
the real background of the theoretical model. We can see 
that there is a remarkable agreement between the two 
in the whole range of invariant masses. We might con- 
clude from there that the mixed event method is really 
good to reproduce the background. Yet, let us investigate 
with more detail how this shape has been produced. Upon 
renormalization of the generated mixed event distribution 



we reproduce the real background, but we know that one 
is using for sure information not included in the spectrum 
of Fig. ED 

Another way to proceed is to select two independent 
events from Fig. [TJ meaning that each one separately ful- 
fills \p n a + p 7 | < 500MeV/c, and then reconstruct the 
invariant mass of Eq. (|24p for the mixed events, impos- 
ing also the cut of Eq. (|30| to the pair of events of the 
mixing. We call that method II. This would correspond to 
a ME reconstruction from experimental events that have 
been filtered with the \p^o +p 7 | < 500 MeV/c condition, 
which imposes a certain correlation, which might be unde- 
sired, in the events chosen for the mixing. The results can 
be seen in Fig. [3J Now we normalize the background at 
low invariant masses where it is maximum. Then we ob- 
serve that in the rest of the invariant mass region there is 
a clear disagreement of this new mixed events background 
with the real one. Hence, the result has been a very poor 
reproduction of the real background by the mixed event 
method II. 

Going back to method I, and in order to understand 
what is really happening, we have conducted another test. 
Let us realize that the mass distribution of Fig. [TJ is expo- 
nential and there are three orders of magnitude difference 
between the strength of da/dM^a^ at low and large invari- 
ant masses. From pure statistics it looks quite logical that 
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Mixed events: Method I 
Mixed events: Method II 
co(782)-signal & real background 
Exp. data (CBELS A/TAPS) 




500 600 700 800 900 1000 1100 

M o [MeV] 

k y 

Fig. 3. The effect of the momentum cut on the mixed events before and after the mixing (method II). 



if we take two independent events to reconstruct the mixed 
event ir°j invariant mass, these two events belong to the 
region of the spectrum that has larger cross section, even 
if the mass that we obtain corresponds to the large invari- 
ant mass where da/dM^o^ is small. In other words, it is 
perfectly acceptable that the background that one obtains 
in the large invariant mass region is largely determined by 
events far away from this region, sitting at much lower 
individual 7r°7 invariant masses. If this were the case one 
would be attributing the background in the high invariant 
mass region to different reactions than those responsible 
for it and hence one would be distorting the physics of 
the process. Certainly in such a case there would be an 
interesting side effect: the distribution obtained with the 
mixed event method at large invariant masses would be 
largely insensitive to the actual background contributing 
in that region. In this case the ME method would thus 
render a background in this region that has nothing to do 
with the actual one. 

In order to illustrate more dramatically the problem, 
we change arbitrarily the background of our model at large 
M„.o 7 by imposing 



where /(M 7r o 7 ) is a distortion factor. We consider two 
sharp cuts with 



(1) /(M w o 7 ) = Q 



1 for M 7r o 7 < 850 MeV 
for M 7r o 7 > 850 MeV 



and 



(2) f(M«o„) 



for M, 
for M, 



7T U 7 



< 750 MeV 
> 750 MeV 



(32) 



(33) 



|T| 2 |T| 2 /(M w o 7 ) 



(31) 



The distortion factor in Eq. (|32j> cuts off the background 
at higher invariant masses beyond the w-signal and (f3"3"f 
removes both the signal and the background. 

In Fig. we show the results for the background from 
the ME method, with method I and the real data, com- 
pared with those obtained with the distorted spectra of 
Eqs. f3"2"|) and (|3"3")) . Note that the sharp cut off becomes in 
the figure a smoother fall down because we implement the 
folding of the 7r°7 invariant mass with the experimental 
resolution of [T3] of 50 MeV. What we see in Fig.0]is that 
the ME method output barely changes with the actual 
background from Eqs. (j32!33[) , or in other words, that the 
ME method is unable to produce the actual background. 
It produces a background largely tied to the events at 
low invariant mass and is unsuited to produce a realistic 
background in the region of large M T o 7 masses. 




M o [MeV] 

k y 

Fig. 4. The mixed event background with Method I for different input invariant mass distributions. 



In order to test the former suggestion that the back- 
ground at large invariant masses in the ME method is tied 
to events at low invariant masses, we construct the cor- 
relation matrix C(M 7T j(ini),M^[fin)) where M ni (ini) 
refers to any of the two events used in the mixed event 
simulation and M V1 (fin) refers to the ME final invariant 
mass determined through Eq. (T24l . 

In Fig. [5] we plot the correlation function. We have 
taken M ff7 (/m) — 800 MeV. Then keeping this variable 
fixed we plot on the y-axis in an arbitrary scale the num- 
ber of events (summing the two events used in the mix- 
ing) that would have a certain M vl [in%). As we can see, 
the initial events used in the mixing accumulate in the re- 
gion where M wry (ini) is about 400-500 MeV. This certainly 
is distorting the physics of the problem, since the back- 
ground associated to the region M vl {fin) = 800 MeV is 
generated after mixing by events around 400 MeV, where 
the origin of the background is quite different from the 
real one around M„^(fin) = 800 MeV. 

6 Summary 

In summary, what we have seen is that due to the pecu- 
liar shape of the background in the present process and the 
fast drop as a function of M„- 7 , the mixed event method 
is unsuited to provide an even qualitative reproduction 



of the real background of the process. Even if a first run 
seemed encouraging because it gave a good reproduction 
of the background, further insight into the method re- 
vealed its flaws since we could prove that different meth- 
ods to do the cuts gave rise to very different mixed events 
background. Further we could prove that the results pro- 
vided by the mixed event method were practically insensi- 
tive to the value of the real background beyond 750 MeV, 
to the point that we could take any arbitrary background 
as input in that region and the mixed event method would 
always provide the same background, with no resemblance 
to the one that it was supposed to reproduce. The study of 
the correlations of events gave us an explanation for this 
finding, since we saw that even at large invariant masses, 
the events of the mixing that generated the final M I7 
were collected for the region of M ff7 (mi) around 400-500 
MeV. In that region the origin of the background is very 
different from the one for the real events at large invari- 
ant masses, such that the mixing event method not only 
produces an unrealistic numerical background, but gets it 
from physical processes quite different from those respon- 
sible for the real background at large invariant masses, 
thus grossly distorting the physics of the process. 
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